Nonequilibrium Transport in Quantum Impurity Models 
(Bethe-Ansatz for open systems) 
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We develop an exact non-perturbative framework to compute steady-state properties of quantum- 
impurities subject to a finite bias. We show that the steady-state physics of these systems is captured 
by nonequilibrium scattering eigenstates which satisfy an appropriate Lippman-Schwinger equation. 
Introducing a generalization of the equilibrium Bethe-Ansatz - the Nonequilibrium Bethe-Ansatz 
(NEB A), we explicitly construct the scattering eigenstates for the Interacting Resonance Level model 
and derive exact, nonperturbative results for the steady-state properties of the system. 
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The recent spectacular progress in nanotechnology has 
made it possible to study quantum impurities out-of- 
equilibrium The impurity is typically realized ex- 
perimentally as a quantum dot, a tiny island of electron 
liquid attached via tunnel junctions to two leads (baths 
or reservoirs) held at different chemical potentials. As a 
result of the potential difference, an electric current flows 
from one lead to another across the quantum impurity. 
The description of such an out-of-equilibrium situation in 
a strongly correlated system is a long standing problem 
and has not been given even in the simplest case of when 
the system is in a steady state. 

In a steady state the system properties do not change 
with time even when out of equilibrium. Such a state is 
reached only under special conditions: each lead needs to 
be a good thermal bath and infinite in size (equivalently, 
the bath level spacing tends to zero.) It then follows that 
particles transferred from one lead to another dissipate 
their extra energy in the lead and equilibrate |2|. 

There are two equivalent ways, time-dependent and 
time- independent, to describe the establishment of a 
steady state in the system. In the time dependent pic- 
ture the quantum impurity is coupled to the two baths 
in the far past, to, and is allowed to evolve adiabati- 
cally under the conditions described above. After a suffi- 
ciently long time, at t = say, a steady state is reached. 
Two elements are required to fully determine the sys- 
tem: a hamiltonian to describe the time evolution and 
an initial condition, po, describing the system in the 
far past. The hamiltonian is chosen to be of the form, 
H{t) = H + e vt Hi, where H describes the two free leads 
(thermal baths), H\ is the interaction term between the 
leads and the quantum impurity, and rj an infinitesimal 
parameter, small enough to ensure adiabaticity yet large 
compared to the level spacing in the leads. The initial 
condition is typically given by, 



Po = 



e -0(H o -j: i tiNi) 
TVe-^ffo-EiMiJVi) 



(1) 



t > t , the system is described by a density ma- 
trix p{t) = T{e iJ t» dt ' H ^}p T{e- iJ 'o dt '^*')}, and the 
properties of the system are calculated in the usual man- 
ner, (0(t)) — Tr{p(t)0}. The establishment of a steady 
state follows, in this language, from the existence of the 
limit to — > — oo with the expectation value becoming 
time-independent, (O) = Tr{p s O} where p s = p(0). 

At T = the description simplifies. The initial con- 
dition is typically given by a particular eigenstate of 
Ho, |3>) baths, describing the baths, each with its own 
chemical potential pi. The steady state is then ob- 
tained by evolving the initial state in time, |4 f ) s = 
T{e*^-«= dt ' H{t,) }\$) baths- The expectation values in the 
steady state are computed from, 
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An equivalent way to describe a non-equilibrium 
steady-state is by means of a time-independent scatter- 
ing formalism. The state \^} s is obtained as an eigen- 
state of the full hamiltonian H = Ho + Hi , satisfying the 
Lippman-Schwinger equation, 
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with pi and iV, the chemical potential and number op- 
erator for particles in lead i. Subsequently, at times 



with |$) baths - the incoming state. The scattering eigen- 
state |^) s can be viewed as consisting of incoming parti- 
cles (the two free Fermi-seas) described by \$)baths and 
reflected outgoing particles given by the second term in 
the above equation. Once again two elements are re- 
quired to fully determine the system: a hamiltonian and 
a boundary condition, \<&) baths, which describes the scat- 
tering state far from the impurity. Note that previously, 
in the time-dependent picture |$) baths played the role of 
an initial condition rather than a boundary condition. 
The finite temperature description in this formalism is 
obtained by summing over scattering states weighted ac- 
cording to the Boltzman weights of the corresponding 
incoming states. 

The construction of such eigenstates is a formidable 
task in general. We shall show, however, that it can 
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be carried out for a class of integrable impurity mod- 
els that includes the Interacting Resonance Level Model 
(IRLM) and the Kondo Model. The Bethe Ansatz so- 
lution of these integrable models in equilibrium has led 
to a full understanding of their thermodynamic proper- 
ties. It is based on solving the hamiltonian of a closed 
system, typically with periodic boundary conditions. We 
shall present in this letter a significant generalization of 
the Bethe Ansatz approach to open systems with bound- 
ary conditions imposed by the leads. This approach, 
the Non-equilibrium Bethe- Ansatz (NEBA), allows us 
to construct the fully interacting multi-particle scatter- 
ing eigenestates and compute non-equilibrium transport 
properties, extending Landauer's original approach pj. 
We remark here that our approach differs significantly 
from the recent interesting work by Konik et al. Q who 
also used integrability to compute transport. In contrast 
to their work we model the leads as free Fermi seas rather 
than coupling the chemical potentials to dressed excita- 
tions. 

We focus on the IRLM at T = and defer treat- 
ment of other models to later publications. The IRLM, 
ffiRL = Ei =1 , 2 ,fe e*^«E + £d<#d + E i= i, 2 ,fc(4 d + 
h.c.) + 2J7Ej =1 2 % p Y-c&i&cftd, describes a resonant 
level, edcftd, coupled to two baths of spinless electrons 
via tunneling junctions with strength t. There is also 
a Coulomb interaction U between the the level and the 
baths. The model is closely related to the anisotropic 
Kondo model with the charge states rid = 0, 1 play- 
ing the role of spin states, and td playing the role of a 
local magnetic field. 

Performing some standard manipulations for impurity 
models: expanding in angular modes around the impu- 
rity, keeping only the s- modes, unfolding the model, and 
linearizing around the two Fermi points we have, 



IRL 



[ dx4{x)d^i(x)+e d dU (4) 
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The model thus obtained is a renormalizable field the- 
ory which requires introduction of a cut-off procedure 
to render it finite. The values of the bare parameters 
U, e d , t will be renormalized as the cut-off is removed to 
yield a physical theory. The renormalized theory cap- 
tures the universal physics - where voltages and temper- 
atures are small compared to the cut-off (bandwidth D). 
The chemical potentials for the leads are not included in 
the Hamiltonian. Instead, they enter as nonequilibrium 
boundary conditions specifying the scattering-state far 
from the impurity. 

We wish to calculate the expectation values in the 
steady state of the dot occupation, fid — d^d, and the 
current operator, / = §Ei=i 2(~l) : ' i /v / 2(V'j(0)d — h.c), 



the latter deduced from I = ±[(iV"i -N 2 ),H]. 

To construct the scattering states we shall use a new 
Bethe- Ansatz technique which, unlike the traditional ap- 
proach based on closed systems and periodic bound- 
ary conditions, allows the determination of a state by 
boundary conditions imposed asymptotically. We shall 
build the many body scattering state using single parti- 
cle scattering states that incorporate the boundary con- 
ditions. It is convenient to introduce the symmetric/anti- 
symmetric basis defined by ijj e / (x) — -^(ipi(x) ± 
"02 (a^)), in terms of which the hamiltonian separates 
into even and odd parts, H e = —i J dx ipl(x)dip e (x) + 
Uipt(0)ipt(0)dU + t($t(0)d + h.c.) + e d dtd, and H a = 
-ij dxtpl(x)dip (x) + Ui>l(0)i/)t(0)<tfd. The boundary 
conditions, however, are imposed in the physical basis, 
"0i/2j requiring appropriate combinations of both the 
even and odd sectors. Both hamiltonians conserve the 
number of particles: H e commutes with N e = N e + Nd 
and H a with N Q . The construction proceeds by consider- 
ing the AT-particle sacttering solutions, N = 1, 2, • • • It is 
important to note that in doing so we have imposed cut- 
offs on the theory. Only upon taking the limit iV, L — > oo 
is the field theory regained and the results become uni- 
versal. 

The single-particle eigenstates of the model take the 
form, 



dx [A{g p {x)^l{x) + e p S) + Bh p (x)^ o (x)}\0) (5) 



with |0) the empty vacuum and A and B arbitrary con- 
stants chosen to satisfy the nonequilibrium boundary 
conditions. We are interested in two solutions, labeled 
±, to the Schrodinger equation for these eigenstates, 
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with e p = tg p (0)/(p— 64) and S p = 2arctan 
writing these states, we have chosen the regularization 
scheme: 8(±x)5(±x) = |<5. Note, that we take h p (x) to 
be discontinuous at zero. This unorthodox choice of so- 
lution is allowed by the linear derivative. Theories with 
linear derivatives - as realized by Dirac long ago - are 
implicitly many-body theories. To calculate physical ob- 
servables, one must first fill the Fermi-sea from a lower 
cut-off D to the Fermi energy. Since the universal many 
body physics is only sensitive to the amplitude of h p (x) 
before and after the impurity, renormalizability implies 
that physical observables are insensitive to the choice of 
discontinuities in h p (x) ||- 

We construct two kinds of single-particle scattering 
states, namely those with incoming particles from lead-1, 
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LEAD 1 SCATTERING STATE 



LEAD 2 SCATTERING STATE 



FIG. 1: There are two types of single particle scattering 
states. In type-1 scattering state, an incoming electron in 
lead 1 is scattered by the impurity and can hop on to either 
lead 1 or lead 2. Notice there is no amplitude for an electron 
to initially be in lead 2. In a type-2 scattering state, the role 
of the two leads is reversed. 



\lp), and those with incoming particles from lead-2, \2p), 
with p the momentum of the incoming particle. Choosing 
A = B in cq(J2Jl, the amplitude for an incoming particle 
from lead-2 vanishes and we get, 



\lp) = / dxe ipx 



T — K -{[28{-x) + {e^ + l)8{ X )]4 
[(e i5 " - l)0(z)]v4) + V2e p d^6(x)] |0) 



Conversely, choosing A = — B, the amplitude for an in- 
coming particle from lead-1 vanishes and we get the state 
\2p), given by the above expression with ip\(x) and V^ 2 -) 
interchanged. It is convenient to introduce the opera- 
tors, a\ /2p (x) = g p {x)il>l(x) ± h^(x)ipl(x) + e p 8{x)d) , 
in terms of which the scattering states are 1 1 / 2p) = 
(x)|0). The single particle scattering eigen- 
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l/2p 



J dxe 

states are depicted in FIG. 1. 

The most general two-particle eigenfunction is of the 
form, 



[Ag(xi, x 2 )i)l(xi)ipl (x 2 ) + Ch(xi, X2)tpt(xi)ipl(x 2 ) 
Bj(x 1 ,x 2 )i>i(x 1 )i>Ux 2 )}\0) 
+ I [Ae[x)i>l{x)d) + Bf(x)4(x)J]\0) (7) 



with, 



2g(xi,x 2 ) 
2h(xi,x 2 ) 
f\xi,x 2 ) 



= 9p{x\)gk{x 2 )Z(xx - x 2 ) - 

= h p {x\)hk{x 2 )Z{xi- x 2 ) - 

= g p {xi)h%{x 2 )Z{xi - x 2 ) 

+ (-l) ab g k (x 1 )h b p (x 2 )Z(x 2 

a/b. 
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with a, b = ±, g p (x) and h p (x) being the single particle 
eignefunctions eq® and Z{x x - x 2 ) = e Mp*)s g n( Xl -x 2 ) 
with e 2 ^^ = (i + % h r\ )/(i " T t/~9 )■ Thc 

V 2 k+p-2e d 'I V* 2 k+p-2e d ' 

constants A, B, and C are determined by the nonequi- 
librium boundary conditions. In this solution, we made 
use of freedom afforded by the linear dispersion to choose 
the two-particle S-matrix between all electrons to be the 
same. This allows us to easily generalize the construction 



to iV particle wave functions yielding the fully-interacting 
scattering state, 



N ± +N 2 



|*) a = I dxV{x 1 ---x N )Y[a{ p Jx u ) J[ a\ Pv {x v )\Q) 

u=l v=N 1 + l 

y s (X! ■■■X N ) = B^i Pi^e^x* v(Ps, P t)sgn(x s -x t ) ^ 

Recall that for \^f) s to describe a non-equilibrium 
steady-state the incoming particles in the region 
{ x j} < must be described by \&) baths- In 
the coventional Fock basis |<6) baths is given by: 

Y\uU elT - uPuXu Uvlthli elT/vPvX ^ with the Fock mo- 
menta {pj} satisfying: — D < p u < p,\ and —D < p v < 

fi 2 . Notice, however, that in ^ s (xi ■ ■ ■ xn) there is a 
two particle S'-matrix, S — e 2tip ( p > k \ between incoming 
particles in each lead though the particles are free elec- 
trons. The presence of this non-trivial S-matrix forces 
a choice of a different, "Bethe-Ansatz" , basis of eigen- 
states for the free Fermi seas in the leads, inherited from 
the interacting model when the coupling to the impurity 
is turned off [2J. In order impose the boundary condi- 
tion in the Bethe-Ansatz basis, the incoming particles 
"Bethe-Ansatz" momenta {pj} in |\I/) S , thus far undeter- 
mined, must be appropriately chosen. This is done below 
by solving "free-field" Bethe-Ansatz equations. 

The steady-state current and dot occupation in the 
non-equilibrium steady-state is computed from eq(|2"|l 
with O the appropriate operator and given by eq(??). 
When computing this expectation value, one must take 
the system size, L, to infinity (recall that no steady state 
can be reached otherwise). In this limit, scattering states 
of type 1 and 2 are orthogonal and we find, 
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(n d )s = E 



^( Pu -e d ) 2 + A 2 v= ^ +i(pv -e d y + A* 
A w i±f ' A 



+ E 



^{Pu-etf + A 2 v ^ +i ( P 



A 2 ' 



with A = t 2 /2. In the non-interacting case, U = 0, 
imposing the boundary condition in the thermodynamic 
limit, the sums are replaced by integrals over pi - product 
of the density of states v = 1/2tt and the Fermi-Dirac 
function - that describe the distribution of momenta in 
each lead, e.g. at T — 0, Pi{p) = -^9(k l — p), with fc* set 
by p 1 . We obtain the standard RL results 



0). 



\ n d) 



dp[pi{p) - p 2 (p)} 



(p - e d ) 2 + A 2 
dp \pi(p) + pa(p)] , _ , 2 + A2 (10) 



In the interacting case, however, pi(p) are no longer 
Fermi-Dirac distributions. As explained above, the pres- 
ence of the non-trivial S'-matrix requires the distribu- 
tion in each lead to be obtained by solving a set a free 
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Bethe-Ansatz equations. For T = 0, the distributions, 
Pi(p),i = 1,2 satisfy, 



Current vs Voltage (D=50, Delta=0.5) 



1 f ° 

pi(p) = ir0{k l o -p)- Y, / ic(p,k) Pj (k)dk 

P2(p) = ^0{k 2 o -p) - J2 t° ^ k) Pj (k) dk 

with k l ,i = 1,2 upper bounds on the distributions of k 
set by the chemical potentials pi (we choose k\ > fe^), 
^b, k) = %{k - e d )/[(p +k- 2~e d f + ^-(p- kf] and 
(-a = — U X)i=i 2 T 1 - The equations need be solved 
in the presence of a cut-off £>, — £) < k. For T > 
one needs to solve the corresponding finite temperature 
Thermodynamic Bethe Anatz equations. 

In the U = oo limit, these equations can be solved by 
a standard, if tedious, Wiener-Hopf method yielding the 
results (I). = A (Zi) 



tan 



l 

2tt 



'IV 
v A , 



tan 



-1 /£]_ 



tan 



-1 M2-£d 



tan 



-1 M2 



and 
where 



T fc = L> (§) "+< with e^^) 



Tk is a new 



low energy scale in the problem, related to the Kondo 
temperature in the anisotropic Kondo model. It is held 
fixed as the cut-off and U are sent to infinity. 

More generally, these equations must be solved numer- 
ically with the bandwidth, D, much larger than all pa- 
rameters in the problem to ensure we are in a universal 
regime. In FIG. 2, we plot our results for the current as 
a function of voltage for various values of U. Notice the 
current is non-monotonic in U with a duality between 
small and large U . This duality holds for all A. We also 
plot the distribution function in lead 1, pi, as a function 
of momentum for various voltages. The strong depen- 
dence on momentum and voltage of these distributions is 
a hallmark of the nonequilibrium physics. 

For pi = P2 the system reverts to equilibrium and 
our construction can be compared with traditional Bethe 
Ansatz approaches. There is no current in this case, 
and the dot occupation can be obtained from the im- 
purity energy, given at zero temperature by Ei mp = 
J dpp(p)S p with p(jp) determined by the TBA equa- 
tion. Hence, ignoring d ed p{p), which is suppressed by 
N/L, we have (n) d = d Cd E imp = J dpp(p)d ed S p . Since 
d e J p = 2A/((p - e d ) 2 + A 2 ), it coincides with ea ((Tfj|) 
when pi(p) = p 2 (p) = p{p). 

In conclusion, we have presented an exact solution of 
a strongly correlated impurity model out of equilibrium. 
The solution is given in terms of the scatteringing states 
that characterize the nonequilibrium steady state. The 
generalization to finite temperature or to more than two 
leads is straightforward. The latter allows the compu- 
tation the nonequilibrium density of states which is 
of experimental interest. We believe the framework we 
introduced is very general and can be applied to most 
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FIG. 2: Here we show the current as a function of voltage for 
various U with p,i/ 2 = £d ± V/2 for fixed bandwidth D and 
A. Note,the current is not monotonic in U. We also show the 
distribution in lead 1 as a function of momentum for various 
voltages where with out loss of generality we take fc* = 0. 



integrable models. Thus far we have constructed current 
carrying scattering states for the Anderson and Kondo 
models, though we do not know a general criterion for 
the framework's applicability. Neither do we have a clas- 
sification of the operators whose non-equilibrium expec- 
tation values can be calculated. 
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